Positronium signature in organic liquid scintillators for neutrino experiments 
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Electron anti-neutrinos are commonly detected in liquid scintillator experiments via inverse beta decay, 
by looking at the coincidence between the reaction products, neutron and positron. Prior to positron 
annihilation, an electron-positron pair may form an orthopositronium (o-Ps) state, with a mean life 
of a few ns. Even if the o-Ps decay is speeded up by spin flip or pick off effects, it may introduce 
distortions in the photon emission time distribution, crucial for position reconstruction and pulse shape 
discrimination algorithms in anti-neutrino experiments. Reversing the problem, the o-Ps induced time 
distortion represents a new signature for tagging anti-neutrinos in liquid scintillator. 

In this paper, we report the results of measurements of the o-Ps formation probability and lifetime, 
for the most used solvents for organic liquid scintillators in neutrino physics (pseudocumene , linear alkyl 
benzene, phenylxylylethane, and dodecane). We characterize also a mixture of pseudocumene +1.5 g/l 
of 2,5-diphenyloxazole, a fluor acting as wavelength shifter. 

In the second part of the paper, we demonstrate that the o-Ps induced distortion of the scintillation 
photon emission time distributions represent an optimal signature for tagging positrons on an event by 
event basis, potentially enhancing the anti-neutrino detection. 

PACS numbers: 13.15.+g 14.60. Pq 29.40. Mc 36.10.Dr 



INTRODUCTION 

Electron anti-neutrinos are produced in f3 decays of natu- 
rally occurring radioactive isotopes in the Earth, represent- 
ing a unique direct probe of our planet's interior. Recent 
results by Borexino [lj and KamLAND [2J have been com- 
bined to provide the first observation of geo-neutrinos at 
5 cr PI- Also nuclear reactors provide intense sources of 
antineutrinos, which come from the decay of neutron-rich 
fragments produced by heavy element fissions. In the past, 
reactor anti-neutrinos allowed to constrain the neutrino os- 
cillation parameters in the Solar sector [4|. Presently, sev- 
eral experiments (Double Chooz |5j, Daya Bay |6J, RENO 
|7j) exploit reactor anti-neutrinos to investigate the last un- 
known neutrino oscillation parameter, 813. Moreover, re- 
actor anti-neutrinos are directly linked to the fissile isotope 
production and consumption processes and can be studied 
for monitoring the nuclear plants |8J. 

One of the main advantages in detecting electron anti- 
neutrinos is the clean signature provided by the inverse beta 
decay reaction: 

v e +p n + e + (1) 

by looking at the neutron-positron coincidence. 

In organic liquid scintillators, neutrons are captured 
mainly on protons, and identified looking at the charac- 
teristic 2.223 MeV 7-ray emitted in the reaction: 

n+p^d + j (2) 

To enhance the neutron detection, scintillators can be 
loaded with high neutron capture cross section materials, 
like Gadolinium. The neutron mean capture time varies 
from few tenths up to hundreds of microseconds. 



Positrons interacting in liquid scintillator may either di- 
rectly annihilate with electrons or form positronium (Ps) |9J, 
a bound state of a positron and an electron with two sub- 
states: 25% of positronium populates the spin singlet state 
(S=0), called para-positronium (p-Ps), and 75% the spin 
triplet state (S=l), called ortho-positronium (o-Ps). p-Ps 
annihilates emitting two 7-rays of 511 keV with a mean life 
of 125 ps in vacuum. o-Ps emits three 7-rays, with a total 
energy equal to twice the electron mass, and a mean life in 
vacuum equal to ^140 ns. 

In matter, however, interactions of o-Ps with the sur- 
rounding medium strongly reduce its lifetime: processes like 
chemical reactions, spin-flip (ortho-para conversion at para- 
magnetic centers), or pick-off annihilation on collision with 
an anti-parallel spin electron, lead to the two body decay 
with lifetimes of a few nanoseconds. The surviving three 
body decay channel is typically reduced to a negligible frac- 
tion. 

If the delay introduced by the positron annihilation life- 
time is of the order of a few nanoseconds, calorimetric scin- 
tillation detectors, like Borexino [H5] and KamLAND [TT] , 
are unable to disentangle the energy deposited by positron 
interactions from that released by annihilation 7-rays. In 
these cases, a delayed 7-ray emission induces a distortion 
in the time distribution of detected photoelectrons (pulse 
shape), with respect to a direct annihilation event. Such a 
distortion can affect algorithms based on the pulse shape, 
like position reconstruction and particle discrimination. 

Distortions in positron pulse shapes affect not only anti- 
neutrino detection via inverse beta decay, but also the re- 
construction of j3 + decay events, like those from cosmo- 
genic 10 C and U C radioisotopes, which represent a crucial 
background for the Solar neutrino signal in Borexino and 
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Experiment Scintillator Fluor Dope 



KamLAND [4J 20% PC 
80% OIL 
Borexino [TOJ PC 
LVD [16] Paraffin 
SNO+ [IZl LAB 
Double Chooz [5] 20% PXE 
80% OIL 
Daya Bay [6] LAB 

RENO [7] LAB 



1.5 g/l PPO 

1.5 g/l PPO 
1.0 g/l PPO 

PPO 0.1% Nd 

3-6 g/l PPO 0.1% Gd 

20 mg/l Bis-MSB 

3 g/l PPO 0.1% Gd 

15 mg/l Bis-MSB 

1-5 g/l PPO 0.1% Gd 

1-2 mg/l Bis-MSB 



Table I: Scintillator composition used in present and future 
underground neutrino experiments. 

KamLAND [12] LE]. 

In this paper, we characterize the o-Ps formation prob- 
ability and lifetime for the most popular choices of scintil- 
lator solvents used by present experiments (table |T|): 1,2,3 
trimetilbenzene or pseudocumene (PC, C9H12), linear alkyl 
benzene (LAB, C18H30), phenylxylylethane (PXE, CigHig) 
and dodecane (OIL, C 12 H 2 6)- 

Wavelength shifters are present in scintillators with very 
low concentrations. However, adding a fluor to scintillators 
can change the o-Ps properties. For this reason, we studied 
a typical scintillator mixture, PC + 1.5 g/l of PPO (2,5- 
diphenyloxazole), to observe the effect of the fluor on the 
o-Ps properties. 

In the second part of the paper we discuss the results of 
Monte Carlo simulations in which we quantify the distortion 
on the pulse shapes introduced by the o-Ps and the effects 
on the positron event reconstruction in large volume liquid 
scintillator detectors. 

Finally, we demonstrate that the o-Ps induced distor- 
tion potentially represents a powerful signature for anti- 
neutrinos. 



EXPERIMENTAL APPARATUS 

The positron annihilation lifetime spectra have been mea- 
sured with a fast coincidence system, formed by two plastic 
scintillators, with a time resolution of 0.28 ns at FWHM. 
22 Na is a (3+ radioisotope, with Q-value = 2.842 MeV, emit- 
ting positrons in association with 1.274 MeV 7-rays (99.94% 
branching ratio). The 22 Na source (0.8 MBq activity) was 
prepared by drying a droplet of 22 Na from a carrier free 
neutral solution between two layers, each composed of two 
Kapton foils (7.5 fim thick, 1 cm radius each single layer). 
Kapton is a polymide compatible with scintillator materials 
and where positrons do not form o-Ps states |14l I15j . 

The Kapton-source sandwich is inserted in a glass vial, 
containing the scintillator sample. The vial is positioned 
between two plastic scintillator detectors (Pilot U) each of 
them coupled with a photomultiplier tube. The vial is in 
direct contact with the faces of the two detectors, in order 
to maximize the solid angle. 
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Figure 1: Color online. Positron annihilation lifetime spectra 
for the PC and PC+1.5 g/l PPO samples. 

The system measures the delay occurring between the 
22 Na decay and the positron annihilation: when a plastic 
scintillator is hit by a 1.274 MeV 7, from 22 Na decay, the 
second detector opens a gate of 43 ns in order to detect one 
of the annihilation 7's. The time resolution of each detector 
is ~ 130 ps. 

Discrimination between start and stop events occurs on 
the basis of the different energies of the detected 7-ray. The 
first detector has a lower threshold at 900 keV to detect the 
prompt signal from the 1.274 MeV 7, while the annihilation 
gamma deposited energy is recorded if it falls in the [350- 
500] keV energy range. The second energy cut has been 
optimized to avoid the back scattering of the 1.274 MeV 
7's. 

A constant fraction discriminator generates a fast timing 
signal whenever a 7-ray with the correct energy is detected. 
A time-to-amplitude converter, enabled by the start signal 
from the discriminator, produces a voltage linearly increas- 
ing with time, which stops at the arrival of the stop signal 
from the other discriminator: the signal at the output of 
the time-to-amplitude converter is proportional to the time 
interval elapsed between generation and annihilation of the 
positron. This signal is finally digitized by an analog-to- 
digital converter with 4096 channels. 

The time interval corresponding to each channel is 10.6 
ps, and it has been calibrated with a 60 Co source. The 
resolution of the apparatus (0.28 ns at FWHM) is measured 
by looking at the fast (182 ps) coincidence between the two 
7's emitted by the 207 Bi decay, through the 5/2~ 207 Pb 
excited state. The detector resolution is modeled with the 
sum of two gaussians with different standard deviations. 

DATA ANALYSIS AND RESULTS 

The 5 samples (PC, PXE, LAB, OIL, PC+PPO) of scin- 
tillators have been degassed with nitrogen before each mea- 
surement. For each sample, we repeated the measurement, 
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Figure 2: Color online. Fit (red line) of the positron annihi- 
lation life time spectrum (black dots) for the PXE sample. 

at the ambient temperature, from 3 to 5 times, to take into 
account possible systematics due to changes of the initial 
conditions (e. g. temperature and time of nitrogen strip- 
ping). 

The collected total statistics for each sample is ~ 5 x 10 6 
events. Measured positron annihilation lifetime distributions 
of PC and PC+PPO samples are shown in figure [TJ 

Each distribution is fitted with a three component model: 



F(t) = X (t > t0) 



J2 ^ ■ e- t/Tk + C | (3) 



, fe=l,2 



n 



where Ai and t\ correspond to the effective amplitude 
and mean life of the annihilation and p-Ps components, 
which can not be disentangled; A2 and T2 are the o-Ps am- 
plitude and mean life, respectively; C is the flat distribution 
due to accidental background; to is the detector time offset 
and x is a step function. The fit function F(t) is convoluted 
with the detector resolution, modeled with the sum of two 



gaussians: 



(4) 



centered in the same value but with different resolutions 
(<7i and 02), and where gi (i = 1, 2) are the relative weights 
of the two components (gi + 92 = 1)- The detector resolu- 
tion is dominated by the gaussian with o\ ~ 110 ps (gi ~ 
0.8), while 02 ~ 160 ps. The data modeling package used 
in this analysis is the RooFit toolkit [18J, embedded in the 
ROOT package [19J, based on MINUIT. All the parameters 
in the model are free in the fit. An example of fit for the 
PXE sample is shown in figure [2] 

All the normalized \ 2 , obtained from the fits, are in the 
[0.85-0.98] range. Statistical errors for both mean lives and 
amplitudes are at the order of per mil level. t\ is almost 



Figure 3: Color online. The data points represent the o-Ps 
fraction, measured with 1 to 3 Kapton layers, inserted in a 
Plexiglas medium. Each layer is 7.5 /im thick and 1 cm radius. 
The line is the result of an exponential fit. 



constant in all the measurements: the mean t\ is centered 
in 365 ps, with a root mean square of 8 ps. 

To estimate the fraction of positrons annihilating in Kap- 
ton, we covered the 22 Na source with 1 to 3 Kapton lay- 
ers, identical to those used for charactering scintillators. 
The Kapton-source sandwiches were inserted in a Plexi- 
glas medium, characterized by an o-Ps mean life of ~ 2 ns. 
The o-Ps fraction, as function of the Kapton layers, is fitted 
with an exponential, as shown in figure[3] We estimate that 
20.6 ± 0.2% of positrons annihilate in the double Kapton 
layer (15 /um thick). 

In order to derive the o-Ps formation probability (/ 2 ) in 
the analyzed liquid scintillator samples, we subtract the an- 
nihilation component in Kapton {Ax)- 



h 



A 



K 



(5) 



To estimate the systematic uncertainties, we evaluated 
first the o-Ps mean life weighted averages of all the mea- 
surements for each sample. We fitted then the distribution 
of all the deviations of each measurement from the corre- 
spondent weighted average with a gaussian, in a likelihood 
approach. The same procedure has been applied also for the 
o-Ps formation probability. The resulting systematic errors 
for the o-Ps fraction and mean life are 0.5% and 0.03 ns, 
respectively, and dominate over the statistical errors. Final 
results are shown in table [M] and in figure [4] 

All the samples analyzed in this paper are characterized 
by an o-Ps mean life about 3 ns and a formation probability 
about 50%. The o-Ps mean life and fraction for PXE and 
LAB represent the extremes in the value ranges: PXE (LAB) 
is the solvent where o-Ps has the lowest (largest) probability 
to be formed and the shortest (longest) life time. 

Kino et al. [20J characterized the o-Ps properties on a PC 
sample saturated with IM2 with a similar experimental setup. 
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Figure 4: Color online. Results of the o-Ps probability forma- 
tion and mean life for the analyzed samples of scintillators. 

The resulting mean life (3.03 ± 0.02 ns) is compatible within 
2 a with the one reported in this paper, while we observe 
a formation probability 5.2% larger. This small discrepancy 
can be attributed to a not identical concentration of oxygen 
|20J in the scintillator sample, or to different condition pa- 
rameters, like temperature (the o-Ps formation probability 
as function of the temperature has been discussed by B. 
Zgardzinska et al. |21J), during the measurements. 



o-Ps INDUCED PULSE SHAPE DISTORTION 

Positrons and the following annihilation gammas, even 
after o-Ps formation, can not be disentangled in large vol- 
ume liquid scintillator detectors, like Double Chooz, Borex- 
ino and KamLAND, by looking at the photon arrival times 
in the photomultiplier tubes. Such experiments are, in fact, 
characterized by time constants, like the fluorescence decay- 



time (table III) and the photomultiplier tube time jitters, 
typically longer than the delay between positron and anni- 
hilation gamma emission. 

However, the o-Ps formation can induce a significant dis- 
tortion in the pulse shape of positrons. In the positron 
detection, scintillator molecules are first excited by positron 
interactions, and then by annihilation 7-rays. If annihilation 
passes through the intermediate o-Ps state, the annihilation 
component is delayed, and the overall photon emission time 
distribution (PETD) results as the sum of the two compo- 
nents, as shown in figure [5] for 0.5 MeV positrons in PC + 
Material / 2 r 2 [ns] 



PXE 
LAB 
PC 
OIL 



0.466 ± 0.005 2.74 ± 0.03 
0.542 ± 0.005 3.08 ± 0.03 
0.485 ± 0.005 2.96 ± 0.03 
0.506 ± 0.005 3.04 ± 0.03 



PC+1.5 g/l PPO 0.512 ± 0.005 3.12 ± 0.03 

Table II: Final results for the o-Ps probability formation and 
mean life for the analyzed samples of scintillators. 



e + + y induced photons 
e + induced photons 
Y induced photons 




30 40 50 60 

Photon emission time [ns] 

Figure 5: Color online. Simulation of the photon emission 
time (black line) induced by 0.5 MeV positrons, all forming 
o-Ps, in the PC+1.5 g/l PPO scintillator. The scintillation 
components, due to positron (blue) and o-Ps decay 7's (red), 
are shown separately. 

1.5 g/l PPO, annihilating after o-Ps formation. The cor- 
respondent PETD deformation, with respect to the direct 
positron annihilation case, is dominant in the first 30 ns. 

The distortion is energy dependent, since the amplitude of 
the first component is proportional to the positron energy, 
while the gamma component has fixed energy. To study 
such dependency, we simulated positrons with energies from 
0.1 to 5.0 MeV. The ratio (R) of the PETD mean values 
between the o-Ps and the direct annihilation cases, is plotted 
in figure [6] It is notable that the deformation of the time 
distribution is larger for lower energy positrons, where R can 
reach up to ~1.6. At higher energies, the distortion reaches 
a plateau, with o-Ps PETD mean values shifted up to 10%, 
depending on the scintillator composition. 

In neutrino experiments, the o-Ps induced PETD defor- 
mation can have an impact on several algorithms for the 
event reconstruction and discrimination. For instance, the 
position reconstruction algorithms (e.g. 0. Smirnov [22] ) 
strongly depends on the first (t<30 ns) detected photo- 
electrons. A bias in the position reconstruction can induce 
a further bias in the event energy reconstruction, since the 
light collection on the photomultiplier tubes depends on the 
reconstructed position. Also algorithms for particle discrim- 
ination in liquid scintillators can be strongly affected, since 
they rely on the dependence of fast and slow portions of 
the scintillation pulse on the energy loss of the interacting 

Scintillator n r 2 r 3 Ni N 2 N 3 
[ns] [ns] [ns] % % % 

PC + 1.5 g/l PPO 3.57 17.61 59.9 89.5 6.3 4.2 
PXE + 1.0 g/l PPO 3.16 7.7 34 84.0 12.0 2.9 
LAB + 1.0 g/l PPO 7.46 22.3 115 75.9 21.0 3.1 

Table III: Scintillator decay time constants (ri) and ampli- 
tudes (N 4 ) for P particles for PC+1.5g/l PPO [23], PXE + 
1.0 g/l PPO [25], LAB + 1.0 g/l PPO [25] 
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Figure 6: Color online. Ratio between the mean values of the 
o-Ps and direct annihilation photon emission time distribu- 
tions as function of the positron energy for PXE + 1.0 g/1 
PPO, LAB + 1.0 g/1 PPO, and PC + 1.5 g/1 PPO, assuming 
the scintillation decay constants from table [TTT| 

particle [22 [24] . 



POSITRONIUM SIGNATURE 

Beyond the necessity to minimize systematics in the event 
reconstruction, the o-Ps induced pulse shape distortion po- 
tentially provides a signature for tagging positrons, and 
hence for enhancing the neutrino detection. Borexino and 
KamLAND could directly benefit by a positron/electron dis- 
crimination by enhancing the tagging and rejection of cos- 
mogenic n C /3 + decay, the main background component in 
the pep and CNO solar neutrino energy window [T2] [13]. 
Also and geo-neutrino experiments could be advantaged 
by discriminating cosmogenic 9 Li and 8 He /3-neutron decays 
which mimic the anti-neutrino signal. 

In order to quantify the o-Ps and the direct annihilation 
positron event discrimination, we simulate an ideal spheri- 
cal detector (4 m radius), with 2000 photomultiplier tubes. 
The simulation is Geant4 based. The cathode is a segment 
of a 20 cm radius sphere. The scintillator is PC+1.5 g/1 
PPO, and its properties are taken from F. Elisei et al. [26] , 
The scintillation photon yield is 10,000 photons/MeV. The 
simulation includes optical effects, like Rayleigh scattering, 
absorption and re-emission, and reflections on cathodes and 
on the spherical steel structure. 

We simulated samples of 10,000 positron events, in the 
detector center, directly annihilating or annihilating follow- 
ing the o-Ps formation. Positrons, directly annihilating, 
have the identical pulse shapes of equivalent electrons with 
same energies plus 1.022 MeV, from annihilation gammas. 
Positrons, in fact, thermalize in 300 ps, much faster than 
the characteristic scintillation decay times, and gammas re- 
lease energy mainly through Compton electron cascades. In 
the simulation, positron energy varies from 0.1 to 5 MeV. 
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Figure 7: Color online. Photoelectron time distributions from 
an ideal 4 m radius detector, filled with PC+1.5 g/1 PPO, for 
0.5 MeV positrons directly annihilating (blue) and forming 
o-Ps (red). The structure at ~60 ns is due to optical re- 
flections on detector materials. Both the distributions are 
normalized to 1 event. 

To take into account the photomultiplier tube jitter, we 
smeared the photoelectron time distribution by 1.4 ns and 
we assume an ideal fast electronic chain based on 1 GHz 
Flash-ADC. The photoelectron time distributions (figure[7]) 
are relative to the first detected photoelectron. 

The pulse shapes differ particularly in the peak positions 
and widths. We define a discrimination variable, p, as the 
ratio between the integrals of the distributions between [0- 
18] ns and [18-60] ns. An example of p distribution for 0.5 
MeV positrons is shown in figure|8] Even with a not sophis- 
ticated estimator, the separation between the two samples 
is clear: p from direct annihilation positrons is gaussian 
distributed, centered in 0.76, while positrons following the 
o-Ps decay have a broad distribution of p, depending on the 
o-Ps life time. Because of the already mentioned energy de- 
pendence, we optimized the p < po{E) cut by varying the 
Po(E) threshold with the requirement of 1% direct annihi- 
lation contamination fraction. Results are shown in figure 
[9] The lower efficiency detection at 100 keV, with respect 
to the one at 500 keV, is due to the detector resolution. 

With this technique, we demonstrate that an event by 
event tagging of positrons is possible, with efficiency as high 
as ^25% at 0.5 MeV, equal to the detection efficiency mul- 
tiplied by the o-Ps formation probability. More sophisticated 
algorithms can improve the positronium detection efficiency, 
providing a new signature for anti-neutrino experiments and 
for background rejections in Solar neutrino experiments 



CONCLUSIONS 

In this paper we measured the o-Ps formation proba- 
bilities and life times for the most popular organic liquid 
scintillators used in neutrino experiments. All the scintilla- 
tors are characterized by similar o-Ps properties, with ^3 ns 
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Figure 8: Color online, p distributions for 0.5 MeV positrons 
directly annihilating (blue) and forming o-Ps (red). Both the 
distributions are normalized to 1 event. 
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Figure 9: Color online. o-Ps tagging efficiency (red) based 
on the p < po(E) cut, requiring 1% of direct annihilation 
contamination fraction, for an ideal 4 m radius detector filled 
with PC+f.5 g/1 PPO. The e+ tagging efficiency (blue) is 
obtained multiplying the o-Ps tagging efficiency by the o-Ps 
formation probability. 

mean life and ^50% formation probability. The distortion 
of photon emission time distribution induced by the o-Ps 
intermediate state is a source of systematic in neutrino ex- 
periments, since it affects positron event reconstruction and 
pulse shape discrimination algorithms. 

We further demonstrated that the o-Ps induced distortion 



represents a powerful signature for discriminating a fraction 
of positrons from electrons. This technique can be exploited 
for both enhancing the anti-neutrino detection efficiency 
in 9 13 and geo-neutrino experiments, and for rejecting the 
positron background in Solar neutrino experiments. 
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